Quenched Computation of the Complexity of the Sherrington-Kirkpatrick Model. 
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The quenched computation of the complexity in the Sherrington-Kirkpatrick model is presented. 
A modified Full Replica Symmetry Breaking Ansatz is introduced in order to study the complexity 
dependence on the free energy. Such an Ansatz corresponds to require Becchi-Rouet-Stora-Tyutin 
supersymmetry. The complexity computed this way is the Legendre transform of the free energy 
averaged over the quenched disorder. The stability analysis shows that this complexity is inconsistent 
at any free energy level but the equilibirum one. The further problem of building a physically 
well defined solution not invariant under supersymmetry and predicting an extensive number of 
metastable states is also discussed. 

PACS numbers: 75.10.Nr, 05.50.+q, ll.30.Pb 



Introduction 

The frozen phase of mean field spin glass models dis- 
plays, in the thermodynamic limit, a very high number of 
stable and metastable states. Such a feature is the con- 
sequence of the disorder and the frustration characteriz- 
ing spin glasses and causing the onset of many different 
configurations of spins minimizing the thermodynamic 
potential, organized in the phase space in a rather com- 
plicated way. In order to comprehend the structure of the 
landscape of the thermodynamic potential below a given 
critical temperature a very important theoretical tool is 
the complexity, else called, in the framework of structural 
glasses, configurational entropy, i.e. the logarithm of the 
number of metastable states the system can evolve to. 
The organization of thermodynamic states in complex 
systems is fundamental, e.g., in the understanding of the 
dynamic properties. 

A certain class of disordered mean-field spin-glass 
models, including, e.g., the p-spin interaction models, 
provides a proper description for the main features of 
structural glasses. These are models displaying a sta- 
ble one-step Replica Symmetry Breaking (IRSB) frozen 
phase^^^'^ very well representing the glassy phase of 
amorphous systems. In that case, when a glass-forming 
liquid is cooled down below a certain (glass) transition 
temperature, it loses the ability of visiting different states 
(at least on the observation time-scale considered) and 
this is reflected by the loss of entropy. All the other pos- 
sible states, not selected at the moment of the transition 
of the liquid to a glass, are, anyway, still there from a 
statistical point of view and could, in principle, still be 
reached on different quenches from high temperature or 
on much larger time-scales. The complexity counts the 
many equivalent states that could have been chosen at 
the moment of the quench. Depending on the time scale 
of observation, the complexity is thus a purely dynamic 
concept. 

In the mean-field approximation barriers between dif- 
ferent metastable states of finite life-time become infinite 
in the thermodynamic limit, eventually breaking the er- 



godicity. Complexity can, therefore, be considered and 
computed even at the static level. The presence of many 
metastable states can be detected looking at the relax- 
ation, displaying a dynamical transition with diverging 
relaxation time-scales at a given temperature, the dy- 
namical temperature T^. Below such transition yet an- 
other one can take place, at the static temperature T^, 
this time a proper thermodynamic phase transition to a 
spin-glass phase. ^ The complexity displayed by systems 
whose statics is properly represented by a IRSB solution 
turns out to be extensive. 

For what concerns the Sherrington-Kirkpatrick (SK) 
model, ^ though, that is the mean-field prototype of amor- 
phous magnets (i.e. actual spin glasses), no dynamical 
transition precedes the static one and the frozen phase 
is a Full Replica Symmetry Breaking (FRSB) spin-glass 
phase. Even though the dynamic behavior of systems 
whose static phase is described by a FRSB or a IRSB 
phase seems to be quite different'* the existence of an 
extensive complexity for the SK model has been put for- 
ward since the very early stages of theoretical modeling 
of such materials. 

In the last twenty years basically two different ap- 
proaches have been presented for the behavior of the 
complexity in the SK model. The first one was origi- 
nally introduced by Bray and Moore® (BM); there the 
complexity was analyzed both in the annealed approxi- 
mation, i.e. as the logarithm of the disordered average of 
the number of states, and as quenched average, i.e. the 
average of the logarithm. The second one was initially 
proposed by Parisi and Potters^" who showed that the 
complexity could be obtained by calculating the parti- 
tion function of m distinct real replicas of the system^^ 
and provided the connection with the previous BM for- 
malism by means of a generalization of the two-group 
Ansatz}"^ In that framework they computed an annealed 
solution using the so-called 'unbroken' two-group Ansatz. 
As we will see, such an Ansatz turns out to be equivalent 
to impose a certain supersymmetry on the complexity. 
Indeed, in Refs. [13,14], explicit computation has led to 
a supersymmetric complexity different from the BM one. 
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Over the years it has become more and more evident 
that an important role in the study of the complexity 
of disordered systems is played by the so-called Becchi- 
Rouet-Stora-Tyutin supersymmetry (BRST-SUSY), first 
discovered in the framework of the quantization of gauge 
theories. "'^^ In the context of stochastic field equations it 
can be shown that the generating functional of correla- 
tion functions averaged over disorder leads to an action 
displaying ERST SUSY.^^ ^his formally coincides with 
the average over the quenched random couplings of the 
number of solutions of the mean-field equations for the 
SK model. In such a model, the property of BRST-SUSY 
has been analyzed in Refs. [14,17] and a possible expla- 
nation of their physical meaning has been presented by 
the authors in Ref. [18]. 

In particular, in Ref. [14] evidence for a BRST-SUSY 
complexity functional was brought about. The differ- 
ence with the BM complexity turns out, actually, not to 
be at the functional level, but simply the new solution 
is a second saddle point of the same functional."'^® This 
means that there is just one BRST-SUSY functional (the 
BM complexity) displaying (at least) two saddle points, 
one satisfying the BRST-SUSY^^ and the other one not^. 
Both solutions identify an extensive complexity, comput- 
ing the number of solutions of the Thouless-Anderson- 
Palmer (TAP) equations-'^ in the annealed approxima- 
tion. In Ref. [18] the authors analyzed the two solutions 
and their problems reaching the conclusion that both of 
them had to be discarded as candidate quenched com- 
plexities of the SK model. So far for the annealed level. 

In the present paper we go further analyzing the 
quenched complexity of the SK model in the FRSB 
scheme, modified in order to generalize the approach 
of Alonasson^^ (applicable to systems with IRSB stable 
static phase). The aim is to see whether or not it is possi- 
ble to build a stable (or marginally stable) solution with 
an extensive complexity, such as in IRSB-stable systems 
like p-spin models. 

The program we apply amounts to compute the 
quenched complexity in the FRSB scheme with proper as- 
sumptions. Our computation scheme is formally a partic- 
ular case of the one starting from the three-dimensional 
antiparabolic equation proposed by Bray and Moore in 
Ref. [22]. The solution is, however, more generic than 
the one obtained, under important constraints, by Bray, 
Moore and Young. Indeed, our assumptions are less 
strict then those of Ref. [23] , thus allowing for the con- 
struction of a complexity as a function of the free energy. 
Both that solution and the present one are BRST-SUSY 
invariant. The main results of the present paper were 
anticipated in Ref. 24. 

In Sees. I-II we recall the replica computation for the 
SK model making the Ansatz of infinite replica sym- 
metry breaking in the variational form of Sommcrs and 
Dupont.^^ The SK model is stable only in the limit of 
infinite breaking of the replica symmetry. Even then 
the stability with respect to overlap fluctuations is only 
marginal. In Sec. Ill we formally generalize the Legendre 



transform approach -'^-^ to the case of FRSB stable systems 
and we modify the replica computation scheme accord- 
ingly, adopting a modified FRSB scheme in the spirit 
of Monasson investigation of hidden states for glass-like 
models. In Sec. IV we show the results for the complex- 
ity and the thermodynamic observables obtained both 
by direct numerical solution and by analytical expansion 
around the critical point up to very high order, within 
the modified FRSB scheme. 

The stability of the solution is discussed in Sec. V, 
where we show that, in the SK model, for any value of the 
complexity different from zero, even marginal stability 
is always violated. This means that, requiring BRST- 
SUSY, no exponential growth of the number of states 
with the size of the system at any value of the free energy 
can be obtained. 

What about a BRST-SUSY- ftreafcm^ FRSB quenched 
complexity? Since any supersymmetric solution would 
suffer of the same instability (see Sec. V), in order to 
give account for an extensive complexity this remains 
the only alternative. In Sec. VI we give the most gen- 
eral expressions of the FRSB complexity and of the dif- 
ferential FRSB equations, without imposing any sym- 
metry. Such equations arc, basically, those already pre- 
sented in Refs. [22,23]. We rederive everything using the 
two-group Ansatz formalism, where the breaking of the 
BRST-SUSY is more easily carried out and the BRST- 
SUSY invariant complexity is described by the unbroken 
two-group. From a purely formal point of view, the 
equations are equal to those treated in the study of tem- 
perature chaos by means of coupled real replicas, with 
different boundary conditions. Their resolution is cur- 
rently under investigation. 



THE FRSB SK COMPLEXITY 



The Hamiltonian of the SK model is^ 

1,N 

W = - ^ Jijaiaj 

i<j 



(1) 



where N is the number of Ising spins ai ± 1 and 
are quenched Gaussian random variables with variance 
Jfj ~ 1/N and moan Jij = 0. The over-bar denotes the 
average over disorder. 

Using the replica trick the average free energy can 
be expressed in terms of the spin-overlap replica matrix 
{a,b = 1, . . . , n; n being the number of replicas): 



/3/-P = -^ + ^limiV92^ 

4 4 n^o n ^ 



(2) 



1 . „ / 



lim - log Tr<^ exp I ^ qab ct" 



where a°' is the spin variable relative to the replica a and 
the trace is taken in the replica space. 



3 



In the limit of FRSB the matrix qat is parametrized 
by a continuous non decreasing function q{x), < a; < 1, 
and the free energy functional (2) becomes 



where 



Po{y) = 



4 

+ 00 



l-2g(l)+ / dx q^x) 
Jo 



(3) 



dyPoiy) <l>{0,y) 



V^nqW) 



cxp 



2q{0) 



(4) 



and (/)(0, y) is the solution evaluated at a; = of the Parisi 
antiparabolic differential equation: 



{x,y) = - 



q{x) 



[<P"{x,y)+x P<l>'{x,yy^ 



with the boundary condition 

(t>{l,y) = Tlog(2cosh/3y) 



(5) 



(6) 



We have used the standard notation, denoting derivatives 
with respect to x by a dot and derivatives with respect 
to y by a prime. 



II. VARIATIONAL METHOD 

Functional (3) has to be extremized with respect to 
^(.t). In order to implement the extrcmization one can 
include the Parisi equation (5) and the boundary condi- 
tion at a; = 1, Eq. (6) into the free energy functional via 
Lagrange multipliers P{x,y):^^ 

/ + OD 
dy P{l,y)[<l>{l,y) (7) 
-OO 

-Tlog(2cosh/3y)] - P [ dx f ly P{x,y){^{x,y) 

Jo J -OO 

+ '^-^[4>"{x,y) + x 13 {<t^'{x,y)f]]. 

The functional /^°p is stationary with respect to varia- 
tions of P{x,y), P{l,y), (j){x,y), 4'{0,y) and of the over- 
lap function q{x). Variations with respect to P{x, y) and 
P(l, y) give back equation (5) and its boundary condition 
(6). Stationarity with respect to variations of (/>(a;, y) and 
(^(0, y) leads to a partial differential equation for P{x, y) 
and to its boundary condition at a; = 0: 

P{x, y) = ^ P"i^, y)-x(3 q{x) [P{x, y) ct>\x, y)] ' , 

(8) 

P(0,y) = Po(2/) (9) 
Eventually, variation with respect to q{x) leads to 

/OO 
dyP{x,y)cj>'{x,yf . (10) 
-OO 



The Lagrange multiplier P(x, y) can be interpreted as 
a probability distribution.^^ The value q{x) of the overlap 
function can be associated with the time scale such 
that for times of order states with an overlap equal 
to q{x) or greater can be reached by the system. In this 
picture, P{x,y) becomes the probability distribution of 
frozen local fields y at the time scale labeled by a;.^^ 

Starting from Eq. (10) and from the interpretation of 
P{x,y) as distribution of local fields we notice that the 
function m{x,y) = (j)'{x,y) represents the local magneti- 
zation over the time-scale x. It obeys the equation 

m{x, y) = — Y-m"{x, y) - (3 x q{x) m{x, y) m'{x, y) , 

(11) 



with initial condition 



m{l,y) = tanh(/Jt/). 



(12) 



In terms of the local magnetization one can express 
both the equilibrium magnetization and the field-cooled 
susceptibility of the system: 



eq 



/OO 
dy P{Q,y) m{Q,y) (13) 
-OO 

/OO 
dy P{Q,y)m'{Q,y) (14) 
-OO 

A. Marginal Stability of the FRSB Solution 

In the SK model, with FRSB stable low temperature 
phase, the stability condition can be written as 

/OO 
dy P{l,y) m\l,y) (15) 
-OO 

/OO 
dyP{l,y) [m'{l,y)]HeqT^ 
-OO 



Deriving equation (10) with respect to x we obtain the 
fundamental relation:^^'^^ 



f 

O — ( 



dy P{x,y) [m'{x,y)Y 



1. 



(16) 



and inserting Eq. (16) in Eq. (15) the latter turns out 
to be marginally satisfied. 

Notice that Eq. (16) can be used as well to sim- 
plify the expression of some observables. If the exter- 
nal magnetic field is zero, the minimum value of q, q{0), 
is zero [m(0,0) = 0]. That causes Eq. (9) to become 
P{0,y) = 6{y). As a consequence, Eq. (16) yields 
m'(0,0)^ = 1, m'(0, 0) = ±1, thus implying rUeq = 
and xfc = 1 in the frozen spin-glass phase. 
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III. THE QUENCHED COMPUTATION 

A. Varying the breaking point of the order 
parameter q{x) 

The solution of the self-consistency equation (10) of 
the FRSB solution in the SK model leads to a mono- 
tonic continuous order parameter function q{x),^'^ with 
X G [0, 1] representing a continuous "replica symmetry 
breaking index" or, in the dynamical interpretation of 
the replica computation,*'^^ a "time-scale label" . In the 
absence of external magnetic field q{x) grows from zero 
at X = to its maximum value for x — Xc- From there 
on it remains constant: q{x) = q{Xc) = qiX), Va; G [xc, 1]- 
The Edwards- Anderson parameter qEA = 9(1))^^ also 
called self-overlap, is a physical quantity representing the 
overlap of two configurations belonging to the same spin- 
glass state. '^^ This scenario naturally emerges solving the 
Parisi equation (5) with boundary condition (6) imposed 
at X = 1. The value of is cl decreasing function of the 
temperature, as shown in e.g. Ref. [34]. 

If, instead than at x = 1, we set the boundary con- 
dition of the Parisi equation at a certain Xh < Xc, the 
overlap function extremizing the free energy functional 
develops a discontinuity = q(l\xb) — q{x^) at xi,, 
where q{x'^) is the value of the overlap at the left side 
of Xb while q{l\xb) is the value at the right side (see e.g. 
Fig. 1). In general, the value of q{l\xb) will be different 
from (larger than) qEA- If the boundary condition (6) or, 
equivalently (12), is imposed at any Xb less than one, but 
larger than Xc, no change in behavior occurs: the value 
qEA of the Edwards- Anderson self-overlap is reached any- 
way and the point at which this happens does not shift. 
The above constriiction is a step forward with respect to 
what was done in Refs. [10,23] where the formal con- 
struction was performed keeping such a difference but 
the solutions were computed setting it to zero. 

Such a generalization can be directly implemented in 
the FRSB scheme by imposing condition (18) at x = Xb 
and, at the same time by allowing for a discontinuity Ag. 
This leads to the following free energy functional: 



= - Y [1 - 2a(ik6)] - T (^^^ 

/oo 
dyPo{y) <i>{Q,y) 
-OO 

/oo 
dy P{xb, y) [<t){xb, y) - <l)i{xb, y)] 
-oo 

/■Xf, /-oo 

-/3 / dx j dy P{x,y)\^{x,y) 

Jo J -oo ^ 

'<t>"{x,y)+^x {cj>\x,y)f]] 



-I 



Mxb,y) = 

T>z = dz e 
Pixb,y,z) 



1 



13 Xb 



log / Vz p{xb,y,z) (18) 

(19) 
(20) 
(21) 



2 cosh [py + (3zy/Xqj 
Aq = q{l\xb) - g(Xf^) 



If Xb > Xc, Aq is zero and the boundary condition (18) 
becomes 

4'i{xb,y) = ^log2cosh.Py = (l){l,y) ; V Xb > Xc (22) 

This way we recover the usual unconstrained FRSB so- 
lution of the SK model. 

It is worth noticing that solving Eq. (5) simply with 
boundary condition (6) imposed at at Xb < Xc just leads 
to same result, i.e. to the same discontinuous g(x). The 
choice we make here of introducing Aq already in the 
boundary condition just makes practical computation 
easier. 

For numerical purposes the study of the equation (11) 
for the local magnetization m{x,y) = <j)'{x,y) is better 
suited. For x < Xb, the stationarity Eqs. (11), (12) are 
replaced by 



m{xb,y) = mi{xb,y) = (tanh {j3y + I3z^/Aq^ 

where the average is defined as 

jVz o{y,z) p{xb,y,z) 



jVz p{xb,y,z) 



(23) 



(24) 



q(x) 



0.4 



0.3 



0.2 



0.1 




0.6 



FIG. 1: Overlap functions q{x) at T = 0.6 evaluated resum- 
ming the series expansion in Ta — T (Ts = 1). The dashed 
line represents the usual FRSB q{x) whereas the continu- 
ous Une is q{x) computed with am imposed breaking point 
Xb = 0.3 < Xc = 0.530229(1). In the latter case q{x) reaches 
a plateau at a given xq. Both functions verify the relationship 
q{x)dx =1-T. 



with the following definitions 
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FIG. 2: Order parameter of the FRSB solution obtained by 
numerically solving the Eqs. (8), (11), (10), (26), with modified 
boundary condition for the Parisi equation, imposed at x = xt 
instead than at a; = 1. T = 0.2. For Xb < Xc ^ 0.54173(1) 
a discontinuity occurs. Notice that for xt — 0.0841 (a::iRSB = 
0.0845(2)) the continuous part of q{x) goes to zero in the 
whole interval [0, xirsb]. At Xb = xirsb the free energy of 
the IRSB solution becomes lower then the relative one in the 
FRSB solution. All overlap functions satisfy 1 — J dx q{x) = 



and the self-consistency equation (10) for q{x) is replaced 

by 



q{x) = J dy P{xb, y) (tanh^ + /3z ) (25) 

= q{l\xb) ii X > Xb 

/■°° 2 
q{x) = dy P{x,y) [m(x,y)] \ix <Xb (26) 

J — oo 



cas) of the system interacting through an infinitesimal 
coupling. 

If Zj is the partition function for a given realization of 
the couplings , the starting point for such a computa- 
tion is the replicated partition function Zy"", where n is 
the number of replicas of a single copy of the system and 
m are the real replicas of the system interacting between 
each other through an infinitesimal coupling. Afterward 
we send n to zero, keeping m finite^^'^^ (see also Sec. 
VI). 

Because of the interaction between the m copies the 
resulting solution in the space of n x m replicas is not 
just the one we would get for n replicas repeated m times 
In other words, the values of the elements of the overlap- 
matrix of each copy are not the same one would get in 
absence of the coupling. 

Using the permutation symmetry of the m replicas it is 
possible to perform proper pcrmTitations of lines and rows 
of qab in order to get the highest overlap values grouped 
together in a diagonal block of dimension to. Sending 
n ^ leads to a finite m e [0, 1], representing a sort of 
density of the number of copies. 

Sotting a value Xb < Xc at which we impose the bound- 
ary condition of the FRSB equations, is equivalent to 
cut the refinement of replica symmetry breaking inside 
the diagonal blocks of the Parisi matrix at q{xb). Be- 
yond this point no further breaking occurs and all the el- 
ements of the matrix inside the diagonal blocks take the 
same value q{l\xb) > q{xb). Such a quantity is always 
something less than the Edwards-Anderson self-overlap 
[(/(llxft) < g(l|xc)]. This cane be easily understood since 
it is a sort of average of the monotonously increasing 
overlap function over the q values labeled by a; > a;^. 

Identifying Xb with to, when the 'copies density' to = 
Xb is chosen higher than the critical value x^ the refine- 
ment is complete and the solution coincides with the one 
without discontinuity. 



We notice that Eq. (23) for Xb > Xc, i.e. Aq = 0, 
becomes Eq. (12). 

In Fig. 1 we plot q{x) at T = 0.6 and Xb < Xc by 
means of series expansions around the critical point and 
subsequent Padc resummation. In Fig. 2 we show q{x) 
at T = 0.2 for different values of Xb < Xc, obtained by 
numerical resolution of the Parisi equation. Notice that 
for Xb less than a particular value (we can call it xirsb) 
q{x) becomes a step function. This happens because at 
Xb < a:iRSB the free energy of the IRSB solution becomes 
lower then the relative one for the FRSB solution. All 
overlap functions do satisfy the relation 1 — J dx q(x) = 
T. 



B. The Coupled Real Replicas Approach^"'^^ 
Generalized to FRSB-stable Spin-Glasses 

The equations discussed in the previous section can be 
obtained by considering to equivalent copies (real repli- 



C. Complexity as Legendre Transform of the 
FRSB Replica Thermodynamic Potential 

Making the identification to = Xb and looking at the 
analogy with systems displaying a stable IRSB frozen 
phase we can generalize the approach or Ref. [11] to 
FRSB systems, defining also in this case the complexity 
S as the Legendre Transform of replicated free energy 
P'^'P as 



S(/) = max[/3cc6/-/3x6/-P] 



(27) 



with the free energy / and (3xb as Legendre conjugated 
variables: 



(3xb 
J 



d £ 



dxb 



(28) 
(29) 
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Exploiting Eq. (29), the complexity can be equivalently 
expressed as 



(3 Xf^ 



rep 



dxh 



(30) 



where the relation Xb{f), or vice- versa f{xb), is yielded, 
e.g., by Eq. (29). If Xf, < .Xc, .f'^P is not stationary 
with respect to the breaking parameter Xb and, therefore, 
E(/) 7^ in some /-interval [/o,/th]- As Xb Xc the 
unconstrained FRSB solution is recovered, thus f{xc) = 
/eq and it also coincides with the lower band-edge fo at 
which the complexity goes to zero. 



The solution we are considering satisfies the BRST 
SUSY.">i^ We win show later on (Sec. VI) how this can 
be simply proved, starting from a general formulation 
involving more than one order parameter and reducing 
to the present one when certain relations among them, 
stating the supersymmetry, are imposed. 

In the following we present our results for the complex- 
ity in the FRSBaq Ansatz obtained both by numerical 
resolution of Parisi equation and by analytical expansion 
around Tg. 



Numerical Resolution of FRSB Equations with 
Pseudo-Spectral Integration 



IV. THE FRSB BRST-SUSY SOLUTION AND 
ITS PROPERTIES. 

The explicit expressions for the free energy and the 
complexity as functions of Xb are obtained from Eqs. (29) 
and (30) and respectively: 



(3xb 



[q\l\xb)-q\a 



(31) 



<t>i{xb,y) - -5— {logp{xb,y, z)) 

pxb 



/oo 
dy P{xb,y) 
-00 

l-2q{l\xb)+ [ dxq^ix) 
Jo 



^{Xb) 





4 ^ 

i^Xb^ 



-<^(0,0) 



[q{l\xbf - q{xl 



(32) 



dy P{xb,y) [Pxb 4>i{xb,y) - {\ogp{xb,y,z))] 
In particular, the free energy can be written as 



f = e-Ts 



(33) 



where e is the internal energy density and s the entropy 
contribution inside a single spin-glass state. 

Thermodynamic relations are valid for any Xb and, 
therefore, the expressions for energy and entropy are 



siP,Xb) = (3 



dp 2 



1 



dx (f'{x) 







(3- 



Q yrep 

~ ^b-^ 

Xi OXb 



(34) 
(35) 



-^[l-q{l\xbt 



P^Xb l\q[l-{l-Xb)q{l\xb)] 



+ 



/OO 
dy P{xb, y) 
-00 



Xb 



{logp{xb,y, z)) -y mi{xb,y) 



where we used the fact that the replica entropy 
P'^df'^P/dP is equal to s + E/xb- Notice that, as Aq 0, 
the expression of the entropy tends to the usual FRSB 
one. 



We performed an accurate numerical study of Eqs. (8), 
(11) with imposed breaking at Xb < Xc, exploiting the 
pseudo-spectral numerical integration method used, e.g., 
in Refs. [34-36]. As a comparison we also performed the 
study of the IRSB solution and its complexity (analyzed 
and discussed in Ref. [14]) for the same values of the 
temperature. 

In figure 3 we show the behavior of the FRSBa^ com- 
plexity as a function of the free energy. We observe that 
in the SK model the complexity (30) is non-zero in a 
given interval / e [fo, /th] As shown in the inset of Fig. 
3, the interval of free energies / in which the complex- 
ity is non zero increases as T decreases. This is also the 
dominion of variation of f'^'^^{f)- 

Above the static transition between paramagnet and 
FRSB spin-glass, at temperature Tg, the complexity is 
zero. Below T,, if we look at the behavior of the repli- 
cated thermodynamic potential f^^^ both as a function of 
the state free energy / and as a function of the breaking 
parameter Xb (sec Figs. 4 and 5) we can see that the ex- 
tremal value of the replicated free energy is always at the 
lowest value of the free energy /: Jq — f[xc) = f^^^ixc)- 
This is the value of vanishing complexity and coincides 
with the equilibrium thermodynamic free energy /cq. 

We define x^ax the value of Xb corresponding to /th 
and to the maximum complexity (see Figs. 5-6): In figure 
5 according to the interpretation mutuated from IRSB- 
stable systems, "'^'^ we can sec that for any .Tmax < x < Xc 
the values of / are increasing as Xb decreases: they 
should be the free energy values of the most significant 
states. Their number is exponentially large and equal to 
exp [N'E{f)]. Beyond Xmax the solution is said 'unphysi- 
cal' {f{x) decreases with x). 

On the contrary, as x decreases from Xc to Xmax, the 
total free energy /''*'? decreases. Indeed /''"p is the free 
energy computed by the replica trick and, thus, it has 
to be maximized with respect to variations of the order 
parameter values. Adopting the usual, unconstrained, 
scheme of computation the maximum value of the SK 
free energy, /"^^p = /eq occurs for Xb > Xc- If we compute 
it at any other value of Xb, f'^'^^ is less than this. In figure 
5 the behaviours of both f^P and / as functions of the 
imposed breaking parameter are shown at T = 0.2. 
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Having analized in detail the behaviours of both the 
complexity and the replicated free energy we summarize 
that (a) the maximum of P°^{f) or, cquivalently, the 
minimum of $(/) = / — TS(/), is given by the extremum 
/o of the /-support (the equilibrium free energy in the 
FRSB stable SK model), for any temperature, whereas 
in models where a Td ^ Tg exists (IRSB stable systems) 
such a stationary point is somewhere in the middle of 
the support, going to /o as temperature decreases^"^ '^^ 
(b) the support of the curve S(/) increases from zero as 
temperature is decreased from the critical point (Fig. 3), 
whereas in IRSB stable systems, below T,, the /-support 
decreases as T is lowered. 

Such apparent contradictions can be solved study- 
ing the Ising p-spin model, which has both a IRSB 
and a FRSB phase. ^ Such an analysis will be presented 
elsewhere. 

A closer analysis of the complexity of the SK model 
at a given temperature as a function of the free energy 
rises yet another, more serious, problem. In figure 7 the 
BRST-SUSY quenched complexity is plotted as a func- 
tion of / together with the BRST annealed complexity 
obtained with a IRSB q^xf'^'^^ at T = 0.2. Convexity 
implies that the annealed average must be larger than 
the quenched one and this must hold at any value of /. 
As we can see from the figures 7-8, the threshold value 
/th in the quenched case (vertical line on the right side 
in Fig. 8) is greater than in the annealed one (left side 
vertical line). This leads to a violation of the convexity 
because in the interval of / values between the annealed 
and the quenched threshold free energies the annealed 
complexity is zero while the quenched one is finite. Also 
their maximum values violate the convexity requirement 
as it is shown in Fig. 8. In Fig. 9 the maximum value of 
E is shown versus temperature. In the limit T ^ we 
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FIG. 4: Behavior of /""(/) and $(/) = / - TE(/). The 
two functions converge to /o as f —> fo (the equilibrium free 
energy), where the complexity goes to zero. 
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FIG. 3: Complexity vs. / at different temperatures. The 
interval [/o, /th] decreases as T is increased toward the critical 
temperature = 1. 



FIG. 5: Replicated free energy f'^ and state free energy / as 
functions of x in the region [xmax, Xc] at temperature T = 0.2. 
The maximum value of (/eq occurs for x = Xc, where it 
coincides with the lowest value of /. 



get, for the maximum value of the complexity, the limit 
Um^ S"^(T) = 0.00778(5) (36) 

larger than the value obtained in the annealed approx- 
imation for the SK model (E™^'=(0) = 0.0073)," corre- 
sponding to a IRSB statics in absence of external fields. 



B. Analytical Resolution with Pade Resummation 
of Series Expansion 

As we have seen (figure 1), q{x) is a continuous 

monotonous function for x < X(,. At some xq < xi, it de- 
velops a plateau such that q{x) = q{xo) for a;o < a; < xj,, 
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FIG. 6: Complexity vs. a; at T = 0.2. Smax = 0.105(1), a;c 
0.54173(1). 



FIG. 8: FRSBa? and IRSB complexities. Detail around the 
threshold values of the TAP-state free energy. 
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FIG. 7: FRSBa9 and IRSB complexities versus /. 



FIG. 9: Complexity E versus T for the FRSBa^ solution. The 
dashed part of the curve is obtained by a fit. 



in particular q{x'j^) = q{xQ). Setting Xb near the break 
point of the Parisi solution, Xc, and defining their differ- 
ence as e = Xc — Xb, we derive the following relationships 
valid at any temperature: 



Aq 


= 0{e) 


(37) 


Xo 


= Xc-0{e) 


(38) 


9(1) 


= 9c(l)-0(e2) 


(39) 


q{xo) 


= 9c(l) - 0(e) 


(40) 


q{x) 


= qc{x) + O(e^) X <xq 


(41) 




= fr^ + 0{e') 


(42) 



where the subscript c stays for the case x > Xc, i.e. the 
unconstrained FRSB solution of the SK model. In par- 
ticular, ^0(1)= qsA and fl^P = /eq. 



From the derivative of Eq. (29) we see that 



,9/ 



rep 



i2 irep 



dXh 



Xb- 



dxl 



= 



(43) 



the state free energy / takes its maximal value as a func- 
tion of Xb] as a consequence, smaller values of Xb have no 
physical meaning. Notice that at Xb = a;max the complex- 
ity S(a;b) is also maximal and T,"^'^^ = S(xmax) represents 
the total complexity at that temperature. 

Decreasing Xb, the value q{xo) of the lowest plateau 
decreases; in particular, at Xb = xirsb we have Aq = 
q{l\xb), that is the FRSB solution free energy becomes 
higher than the one of the IRSB function with qo = and 
qi = q(l\xb)- However, a;iRSB < a^max and, therefore, it 
is always outside the "physical" range of Xb^s. 

Wc have computed the function q{x) and the various 
quantities of interest in a series expansion in powers of e 
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FIG. 10: higher curve: qsA as a function of r = 1 — T; 
lower curve: self-overlap g(l|xmax) of the states with ma^ximal 
complexity as a function of r = 1 — T, note the finite slope at 
r = 0. 



and of the reduced temperature t ~ — T (T^ = 1) at 
very high orders through the methods described in Ref. 
[34]. The various series turn out to be asimptotycs but 
they can be resummcd e.g. through Padc approximants 
yielding very precise results in the whole low temperature 
phase. For the total complexity S™'*^(T) at temperature 
T we have the following series expansion: 



E™^(T) = 1 hO (t^°) (44) 

^ ^ 81 81 1215 10935 ^ ^ ' ^ ' 

Higher orders show that this series is strongly asympto- 
tycs (sec Appendix), but Pade approximants of it give 
good results, in particular we have 



lim E'"'«(T) 

T->0 



: 0.0077(2) 



(45) 



in agreement with the result obtained by numerical eval- 
uation of the differential equations of the preceding sec- 
tion. If we compare the quenched complexity with the 
annealed one,^'^^'^^ we can observe that, at any order in 
T, the convexity property is conserved up to the threshold 
energy value of the annealed solution, f^^ (compare with 
figure 8), while it is violated soon after. At the threshold 
fth the difference between the quenched and the annealed 
complexity is 



i:(/,l) - EM) = + O < (46) 

On the contrary, the difference between the maximum 
value of the two complexities is, 



S (r)-S„ (r) = — - — + ^j^+0{T ) 



54675 



that is positive to the leading order, thus violating con- 
vexity. The above expansions up to very high order is 
reported in Appendix. 

The self-overlap of the states with highest free energy 
is always smaller than the Edwards- Anderson parameter 
Qea, Hce figure 10, in particular near T = 0, the self- 
overlap of the states with maximal complexity is 



g(l|a;„,ax) = l-0.4T + O(T2) 



(48) 



while the Edwards- Anderson parameter has a depen- 
dence. 



V. VIOLATION OF MARGINAL STABILITY 

To solve the above mentioned problems the analysis 
must, then, be deepened. The behaviour of the complex- 
ity in itself not being enough to understand the properties 
of the system. We go, then, further analyzing its stability 
in the FRSBa^ scheme of computation. 

A. The Plefka Criterion 

A necessary condition for the physical relevance of the 
solution is the Plefka criterion,^^'^° stating that the quan- 
tity 



(49) 



where the overbar denotes the average over quenched dis- 
order, must be positive or zero. Such quantity coincides 
with the replicon eigenvalue Ar in the replica computa- 
tion. 

In order to compute xp we need the quantity P{x'l, y). 
To this aim one can consider q as the independent vari- 
able instead of x. The discontinuities in q{x) become 
plateaus ofx{q) and viceversa (see Figure 11). The func- 
tions P{q, y) and m{q, y) satisfy the differential equations 



dqm{q, y) 



--jm"{q,y)- 13 x{q) m{q,y) m'{q,y) (50) 



(47) 



dgP{q, y) = -P"{q, y) - P x{q)[m{q, y)P{q, y)]' (51) 

with initial conditions (in the absence of a magnetic field) 
m{quy) = tanh(/3y) (52) 
P%y) = 5{y) (53) 

where qi = q{l\xb). Notice that the previous equations 
do not depend on the derivatives of x{q)] therefore we can 
safely use them also where x{q) is discontinuous. This 
would not be true if we had used q{x) instead of x{q). 
By solving these equations at fixed x = Xb we obtain: 

j dy P{qi,y) tanh'' % = j dy tanh'' y (54) 

^ J dy' GAq{y-y') P{qb,y') cos\f^{py) 
J dy" GAg{y" - y') cosh^'{py" - dy') 
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where Ga{z) = exp[— z^/(2A)] and qi, = q{x^). By 
a simple change of variables the last expression can be 
rewritten as: 



/ 



dy P{qi,y)teinh'' Py (55) 
dy P{qb, y) (tanh'^ P{zyf^ + y) 



I' 



where the average ((■■■)) is defined in Eq. (24), so that 
the Plefka parameter finally reads: 



l-2qi 



(56) 



+ 



J dy P{qb, y) (tanh* (j3y + (3z^fK'^ 



In order to derive the Plefka parameter xp we have to 
compute disorder averages of the single-site magnetiza- 
tion m*^. They can be expressed in terms of the cavity 
field y which is the field induced by spins on the A'' + 1- 
th spin So, averaged in the absence of the spin sq."*^ Intro- 
ducing the distribution function P{y) of the cavity field, 
averaged with respect to the states and to the disorder, 
we have: 



rn'= = j dy P{y) tanh'^ ^y (57) 

It can be shown^^ that the distribution of the cavity field 
at equilibrium is simply related to the function P{x, y) 
introduced to compute the Parisi solution q{x) ."^^ '^^'^'^ 

We now briefly discuss how the same results can be ob- 
tained using the cavity method. According to the Parisi 
picture the average number of equilibrium states with 
given free energy F inside a cluster F is given by 



dM{F) = eyi^[-x{q)l3FT + x(3F]dF 



(58) 



Where Fr is the free energy of the cluster. The cav- 
ity method derives the expression of the free energy and 
the function q(x) by simply using the condition that the 
previous distribution (58) must reproduce itself when a 
new spin sq is added to a system of N spins. If we as- 
sume that there is a complexity curve S(/) we can apply 
the cavity method to reproduce the recipe of Ref. [11]. 
The only difference with the standard equilibrium treat- 
ment is the starting point, i.e. the average number of 
metastable states with free energies near a given free en- 
ergy F* = Nf*, which now reads: 



dJ\f{F) =exp 



ArS(/*) - 



as 

dl 



F -\ 

/=/• 



F 

f=f' 



If we make the formal identification 



dF 
(59) 

(60) 



the derivation of chapter V of Rcf. [41] can be applied 
straightforwardly. There is a crucial difference, however, 
between the two expressions: the first quantity is 0(1) 
for the states at the equilibrium free energy while the sec- 
ond quantity is 0(e^^) for the metastable states. Corre- 
spondingly, in the equilibrium case the final result coin- 
cides with the true equilibrium free energy, while in the 
latter case the result is the free energy of a system in 
which the states a arc weighted with a non-Boltzmann 
weight proportional to e~'^^^° where fix = d'S{f*)/df is 
the Legendre conjugate of /. A Legendre transform must 
therefore be applied to obtain the complexity, just as in 
Ref. [11]. 

In the context of the cavity method the distribution 
P(i/) of the cavity field y inside a state, required to com- 
pute objects like m*^, can be identified with the function 
P{x, y) at the value of x corresponding to the self-overlap 
of the states in the framework of replica computation (see 
Sec. II). The physical meaning of the functions P{x,y) 
and m{x, y) at any x can be understood and their evolu- 
tion operators rederived. The function P{x, y) describes 
the distribution with respect to the disorder of the av- 
erage cavity field y of a cluster labelled by q{x) and is 
computed taking into account the presence of the spin 
So- The function m{x,y) yields the magnetization, av- 
eraged over the cluster, of the spin sq when the average 
cavity field is y. However, the weights assigned to the 
subclusters to compute the averages are not the Boltz- 
mann weigths. Indeed, the averages inside a cluster F 
are defined as: 



{0)i 



-l3xFc, 



{0)c 



-0xF„ 



(61) 



/=/* 



Note that the previous average coincides with the 
Boltzmann averages only at the level of the configura- 
tions corresponding to a; = 1. Thus the magnetiza- 
tion inside a state is actually the thermodynamic mag- 
netization in presence of an average cavity field y, i.e. 
171(1, y) = tanh(/3?/). The peculiarity of these averages 
is that the subclusters with lower free energy no longer 
dominate the average, and all the subclusters have the 
same (infinitesimal) weight. As a consequence, the quan- 
tity / dy P{x,y) m.{x,yY, the average of the overlap in- 
side a cluster, is simply given by q{x), yielding the vari- 
ational equation (10). Had we weighted the subcluters 
with their Boltzmann weight this would not be true since 
few states would have dominated forcing us to consider 
diagonal contributions proportional to the self-overlap. 
The cavity method thus yields the evolution operators 
for P{x,y) and m{x,y) which coincide with the equa- 
tions (8), (11). 



B. Overlap Discontinuity and Replicon Eigenvalue 

In this section we present a simple argument in order 
to show that the FRSB quenched solution must be dis- 
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continuous at Xb in order to implement the recipe of Ref. 
[11] but, at the same time, this discontinuity causes the 
rephcon eigenvalue A/? = a;p to be negative. 

We start considering the expression of the replicated 
variational free energy as a function of q{x) reported in 
Eq. (18). This expression must be extremized with re- 
spect to all parameters but Xb- Therefore, as noted in 
Ref. [25], at the extremum the total derivative with re- 
spect to Xb coincides with the partial derivative: 



in this paper must verify the condition 



(62) 



log / Vz p{xb, y, z) - {\ogp{xb, y, z)) 



where p{xb, y, z) is defined is Eq. (20). The above quan- 
tity is zero when Ag = 0, therefore we need a discontinu- 
ity of q{x) a,i X = Xb in order to compute the Legendre 
transform, Eq. (27). 

In the following we show that the presence of the dis- 
continuity leads to a negative replicon. The discontinuity 
corresponds to a plateau in the function x{q) starting at 
qb and ending at gi. The function E(q) defined earlier 
must be equal to zero at qb and gi but not between these 
points, see figure 11. We work under the hipothesys that 
e = Xc — is small and that the corresponding function 
x{q) is similar to the one of the equilibrium Parisi solu- 
tion (e = 0). In particular, we assume that the function 
x{q) displays a region with dqx{q) 7^ on the left of the 
point q^ . This is precisely the behaviour we have encoun- 
tered in the analisys of the SK model in the preceding 
sections. Our analisys will be valid only at small e, which 
corresponds to the region of free energy near the equilib- 
rium free energy. This region however is quite important 
since it is natural to assume that the function E(/), if it 
exists, is a continuous function of / going to zero at the 
equilibrium free energy. Furthermore, the outcome of the 
present analysis is confirmed by numerical evaluation for 
which it is verified in the whole range of Xb values, as we 
will show in the following. 

It is useful to introduce the fimction E{q): 



E{q) =q- dy P{q, y) m{q, yf 



(63) 



The variational equations impose that E{q) — where 
dqx{q) 7^ 0. Instead, if x{q) is constant, e.g. for 
:< 9 < <lf, we must have E{qi) = = E{qf) at the end 
points but not in the rest of the interval. In other words, 
if there is a plateau in x{q) [which corresponds to a di- 
continuity in q{x)] the equation q = J dy P{q, y) iri{q, y)^ 
must be verified only at the starting point and at the end- 
ing point. We mention that this is sufficient to show that 
the FRSB solution with a discontinuity at Xb considered 



/ dx q{x) — 1 — T 
Jo 



(64) 



as the equilibrium solution actually does. 



Xb 




FIG. 11; Function x{q) and equation E{q). The value a::o is 
the one corresponding to the lower end of the g-plateau: qt = 
q{xo) = q{xt). The ma^ximum overlap value is qi = q{l\xi,). 



The derivatives of the function E{q) with respect to 
q can be evaluated using the equations for P{q, y) and 
to((?, y) written above and expressed only in terms of 
derivatives with respect to y.^^''^^ After some manipu- 
lations we have: 



dqE{q) = 1 - j dy P{q, y)[m' {q, y)Y (65) 
dlE{q) = 2(3x{q) j dy P{q.y)[m' {q,y)f (66) 
dy P{q,y)[m'{q,y)'f 



dlEiq) = 2/3 d,x{q) / dy P{q,y)[m' {q,y)f (67) 



QP^x{qf j dy P{q,y)[m\q,yt 
12pxiq) J dy P{q,y) m' (q, y)[m" (q, y)]^ 
dy P{q,y)[m"'{q,y)Y 



The central point of the argument is the observation 
that the replicon eigenvalue can be expressed as the 
derivative dqE[q) evaluated at qi. Indeed at qi we 
have m = tanh/Jy which satisfies m' — /3(1 — ni?), 
thus dqE{qi) = xp. This property leads immediately 
to the result that each model with a continuous FRSB 
q{x), such as the SK model, is marginally stable, see e.g. 
[25,29] and Sec. II A. 
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The previous expressions allow to study the behaviour 
of the derivatives of E{q) at the point where x(q) is 
discontinuous (see fig. 11). Indeed, since the functions 
P{q, y) and m(g, y) are continuous at qi, we can relate the 
derivatives on the right of q^ to the derivatives on the left 
of qh, which are identically zero. We have, then, 

d,E{q+) = d,E{q^) ^ d,E{q+) = (68) 



The argument can also be generalized to soft-spin sys- 
tems (with soft-spin distribution normalized to (5*^)=!): 
in general at g'l we have 

thus the identification of the replicon with the derivative 
dqE{q), computed a,t q = qi = q{l\xi,), holds once again. 



dlE{q+) = 2p{xo + Ax) j dy P{q,y)[m'{q,y)f (69) 

- j dy P{q,y)[m"{q,y)f 

= dlE{q^) + 2pI^x j dyP{q,y)[m'{q,y)f 
= ciAa; > (70) 
with Ax = Xb — xo- Analogously we have: 

dlE{q+) = -2l3d,x{q^) J dy P{q,y)[m'{q,y)f 

+0{Ax) 

= -C2+0{Ax) < (71) 

The last inequalities in (70) and (71) follow from the fact 
that the function m{q,y) entering J dy P{q,y)[m' {q,y)]^ 
is equal to tanh/3y at first order in Ax and therefore 
at this order the integral is strictly positive. Now the 

relation 

= E{q,) = ^AxAq^ - ^Aq^ + 0{Aq^) (72) 
implies 

Ag = (3ci/c2)Aa; + OiAx^) (73) 

and 

dgEiqi) = ~{Axf^ + OiAx"") < . (74) 

/ C2 

The previous argument can be generalized to any mean 
field spin-glass model with a FRSB phase under the same 
hipothesis for the function q(x) at xi, near Xc, i.e. a 
discontinuity at x = Xf, and a region with dqx{q) ^ 
on the left side of q^ — q{xh). Indeed, in presence of 
multiple p-spin interaction the only change in E{q) is that 
we must substitute the first term q with a function dqf{q) 
which depends on the interaction strength. In general, 
the Plefka parameter xp, or equivalently the replicon A^j, 
can be expressed as dgf{q) — J dy P{q, y) (1 — m^)^, thus 
it can be again identified with dqE{qi). The function 
dqf{q) does not depend on x{q) and therefore it does not 
affect the difference between derivatives on the left and 
on the right side of qi, and the inequalities in Eqs. (70) 
and (71) still hold, leading again to a negative replicon. 



C. BRST-SUSY Subextensive Complexity 

Our most important result is that the quenched FRSB 
BRST-SUSY solution here described must be rejected on 
physical ground at any value of free energy apart from 
the equilibrium one. 

We have checked that the BRST-SUSY solution stud- 
ied in Sec. IV does not satisfy the Plefka criterion [equa- 
tion (49)] , i. e. the replicon eigenvalue is negative as soon 
as Xb < Xc. In figures 12 and 13 we plot the replicon, 
Eq. (56) versus the breaking parameter Xb and the free 
energy /: it turns out to be always negative for Xb < Xc, 
going to zero from below as Xb Xc and remaining zero 
for any value of Xb above the self-consistent point x^- 
Even going beyond the annealed approximation for the 
SK model, ^'''^^ where the relative IRSB static phase was 
not satisfying Plefka's criterion, the condition Ar > 
remains violated. The FRSB BRST-SUSY solution with 
a non zero Aq is not even marginally stable. In order to 
have marginal stability it must be Aq = that means 
Xb > Xc, i.e. the breaking point cannot be imposed but 
it has to be determined self-consistently. 

For such a marginally stable solution, however, S is 
zero, thus there is no way to im,plem,ent the Legendre 
transform approach to compute any physically well de- 
fined complexity. If BRST-SUSY holds, there is no expo- 
nentially large number of states at any energy level above 
/eq in the SK model. 

In terms of free energies, /eq is equivalent to the lower 
band-edge /o = f{xc), since we are computing the stat- 
ics using the right FRSB scheme. This means that, pro- 
vided the BRST-SUSY complexity actually counts some 
extensive number of TAP solutions, they have no phys- 
ical meaning. As discussed in Ref. [18] the violation 
of the Plefka criterion leads also to a mathematical in- 
consistency with a given assumption implicit also in the 
present solution (i.e. the condition _B = 0, sec Refs. 
[9,14,18,22,23]). As we have previously shown in the 
present section, a direct computation indicates that this 
result can be extended to generic mean-field models with 
a continuous FRSB q{x), e.g. the Ising p-spin model be- 
low the Gardner transition (see Ref. [5]). 

Since this solution is unphysical we arc led to the con- 
clusion that, holding SUSY, the number of states of the 
system is subextensive and the states have all the same 
free energy per spin equal to the equilibrium value. Fur- 
thermore all the states verify xp = 0. 
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the quenched complexity 
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FIG. 12: Replicon eigenvalue Ah vs. the breaking parame- 
ter Xb at T = 0.2. Any value of xi, < Xc yields an overlap 
discontinuity Aq causing a negative Ar. 
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FIG. 13: Replicon eigenvalue A^ vs. free energy at T = 0.2. 
The replicon is negative as soon as the free energy landscape 
is probed above its lowest value. 



VI. IS THERE A BRST-SUSY-BREAKING 
SOLUTION? QUENCHED FRSB COMPLEXITY 
IN THE MOST GENERAL SCHEME OF 
COMPUTATION. 



In the previuos part of the paper we have been comput- 
ing the complexity generalizing the recipe of Monasson-*^^ 
for IRSB systems in zero external field to FRSB systems. 
That means that the complexity has been built as Legen- 
dre transorm of the replicated free energy, Eq. (18). In 
other words, the BRST-SUSY has been assumed from the 
beginning. We look now at the possibility of computing 



= ;^ log />«(/) 



(76) 



in a more general formulation, inspired by the Boltzmann 
microscopic definition of entropy. It reduces to the one we 
just employed if particular supersymmetric relations are 
satisfied. The FRSB expression for the quenched com- 
plexity has been presented for the first time, as far as we 
know, in Ref. [22], as 



S, = -A g - u / - A(l - (7) - 



2/32 



(77) 



-I- 



1 

^ dx[p^x)+v*{x)v{x)] +</)(0,0,0) 



where (p{0,0,0) is the "equilibrium", zero local fields 

value of the function <t>{x,hi,h2) solution of the an- 
tiparabolic three-dimensional differential equation 



dhi 



r] 



dhi 

2 



dhidh2 



(78) 

9(6 d(t> 
dhi dh2 



where the derivative with respect to x is indicated by the 
dot. The boundary condition at a; = 1 is 

.^(l,/ii,/i2) =log j dm e^f»"('";'''9'^>^'''(i)>^(i)'''*(i)) 

(79) 

£fuu(TO; u, q, A, A, 7?(1), p(l), r?*(l)) = (80) 
-1 log [27r {P'q - r?(l))] + log + 

[tanh~^ m- {A + p{l)) m + hi]'^ 
2{P^q-V{l)) 

' -|- uf{m; q) + h2m 



+ A- 



/(m; g') = - log 2 -h - log(l - m^) 



(81) 



1 

H — mtanh" m (l — o^) 

The order parameters q, X, A,r]{x),r]* (x), p{x) have 
to be determined self-consistently by extremising (77). 
Since we will reformulate everything in the general- 
ized two-group Ansatz notation we will present the self- 
consistency equations later on, once the linear trasforma- 
tions to that formalism will be applied. We only mention 
that r]{x) is proportional to what was called q{x) up to 
now and that the BRST-SUSY implies a dependency of 
A and A on g and of rj* and p our]. The number q is the 
value of all the elements of the diagonal submatrices of 
the two-group Ansatz, and it is usually set equal to g(l) 
in the unbroken two group scheme. 
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A. BRST relations 

Before going on describing of the two group Ansatz 

wc very shortly recall the relations among order parame- 
ters satisfied by systems that are invariant under BRST- 
SUSY. 

In the BM notation, ^'^^ the macroscopic BRST rela- 
tions, at the level of order parameters, are: 



.18 



-u q 



for the annealed approximation and 



Pab = -^U Qab 

. 2 

nab = — " 9 



(82) 

(83) 

(84) 

(85) 



in the quenched RSB case (see Ref. [44] , where a detailed 
derivation is presented in the notation of Ref. [14]). 



B. Formulation in the two-group Ansatz 

In this section we will first recall the approach of Ref. 
[10] in order to reformulate Eqs. (77 - 80) in the two- 
group formalism. The motivation to perform such a 
transformation is that, from the point of view of the 
BRST-SUSY, imposing the BRST relations is equivalent 
to set equal to zero all order parameters involved but 
q{x). 

The two-group Ansatz was introduced by Bray and 
Moore in Ref. [12] in order to solve the instability prob- 
lem of the replica symmetric solution of the SK model. 
In Ref. [1 1] Monasson showed how the formalism of Leg- 
endre transforms can be applied to mean-field disordered 
models through the pinning of real replicas in a configu- 
ration space extended to m copies of the system. 

Parisi and Potters^" explained how the BM action can 
be obtained through the method of Monasson provided 
that the symmetry between real replicas is broken ac- 
cording to a generalized two-group Ansatz. Being n the 
number of replicas introduced to compute the quenched 
average and m the number of real copies, they analyzed 



1 



lim — log Z 

n— >o n 



mn ^ ext 
/ 



log p,{f)-mpNf 



(86) 



where Z™" is the partition function of n x m copies of 
the system. In terms of the replicated matrix-parameter 
Q, its average is [cfr. Eq. (2)] 



: ext exp 



+ log 




where the indexes a,b = 1, . . . , n while c, = 1, . . . , m. 
The four index matrix can be expressed as the com- 
position of sub-matrices Qab of dimension m x m of 
the form: 



Qab 



Qab 'Sob 
Qab Qab 



The matrices are further parameterized as 



A. 



Qab = Qab ± — ^ + 



y 



Cab 

2y2 



(88) 



(89) 



Furthermore Q^a = 0. In Ref. [10] the last term was 

denoted by Bab/y^- We write Cab ~ "^Bab in order to 
obtain more symmetric expressions in the following and 
to avoid confusion with the parameter B in Eq. (80). 

Evaluating Eq. (86) with this Ansatz. one obtains a 
complexity that can be formally connected to the one 
of Bray, Moore and Young^^ through a given change of 
variables. 

For what concerns the diagonal sub-matrices - cor- 
responding to the annealed case - the transformations 
bringing from the notation of Ref. [22] (BM) to the one 
of the generalized two-group^" (PP) are: 



BM 

q 

A 

A 



PP (two-group) 

Q 



-Q 



y \ C + mA+ —Q 



-m 



(90) 
(91) 

(92) 
(93) 



Writing equations (77)-(81) with the substitutions 
(90)-(93), allows for an immediate connection between 
the breaking of the matrices structure into two groups 
and breaking the BRST-SUSY. Indeed, Eqs. (91) and 
(92) transform into the BRST relations (82)- (83) if we 
set A = C = 0, i.e. if we do not break the matrix struc- 
ture at all {unbroken two-group). On the contrary, setting 
values of A and C different from zero amounts to break 
the BRST-SUSY and leads to independent values of q, 
A and A. 

The changes of variables for the off diagonal terms riab, 
?7*(, and Pab are: 



BM PP (two-group) 

Vab = P^Qab 



Pab 
Vab 



0" Ua6 + 



}ab^ 



(94) 
(95) 



m 



^''[Cab + mAab^'^Qab] (96) 



Setting Aab = Cab = in the two-group construction 

reproduces the generalized BRST relations for the replica 
symmetry breaking case, Eqs. (84)-(85). 
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Is there any FRSB solution breaking the BRST-SUSY? 
If a complexity indeed exists the lower band-edge of its 
dominion in free energy must coincide with the SK equi- 
librium free energy /eq, computed making use of the 
FRSB scheme. Such a value has been found also by 
Bray, Moore and Young: ^'^ starting from the point u = uo 
(such that S(uo) = 0) a BRST-SUSY solution develops. 
Indeed, their solution was based on an assumption equiv- 
alent to the BRST-SUSY and, therefore, it was not the 
quenched analogue of the BM solution (see Sec. VIC). 

As we just underlined , the BRST relations for the or- 
der parameters involved in the computation can be cast 
in a very straightforward way if the notation of general- 
ized two-group Ansatz of Ref. [10] is used. 

In the Full RSB limit the matrices Qab, ^ab and Cab 
tend respectively to the functions q{x), a{x) and c(x) as 
n — > 0. In such limit the expression for the quenched 
complexity (77) (with B = 0) becomes^^ 



>(0,0,0) 

[\x [2 {A^ - a^{x)) -4u{AQ- a{x)q{x)) 
4 Jo 



+2{QC- q{x)c{x)) + (Q^ - (/^(a 



(97) 



where (/)(0, 0, 0) is solution of 



q 
92 



dyl 



dyi 



_dyidy2 dyi dy2 
with boundary condition 



(I){l,yi,y2) = log J dm e^'"" 

= Puf{m: q) - (3m Qyi - ^2) 



(98) 

(99) 
(100) 



+ log 



l-m2 



--log[27r/32(Q-g(l))] 



2 

--m 



C-c{l)-u{A-a{l)) + -{Q-q{l)) 



- 2/?2(g-g(i)) ^'"^""^"^"^ 

-I3^m [(A-a(l))-|(Q-g(l))]}' 

/3/(m; 5) = - log 2 + ^ log(l - m^) (101) 

1 0^ 
-|--mtanh~ m — ^ (l — Q^) 

and where the local fields /ii, /12 in Eq. (78) have been 
changed in yi , 2/2 following the linear transformation: 



yi = — 



/3 



(102) 



2/2 



u hi /12 



Setting 

m = tanh(/3i/i -|- j3z\/ Aq) 
the above formula (99) becomes 

</.(l,yi,y2) = -Y«(l-Q') 



/OO 
Vzp{-u,yi,z) e'- 
-00 



£2(^2, 2/1,2/2) = /3tanh(^t/i + /3^V^g) 

X 



(103) 



(104) 



(105) 



(106) 



Aa 
^^^= +2/2 



+^ tanh(/3yi + ^z^q) {AqAc - (Aa)^) 



where p{w, y, z) is defined in Eq. 20. 

The self-consistency equations for the order parameters 
are: 

q{x) = j dyi dy2 P{x,yi,y2)\dy^(j){x,yi,y2)f (107) 

a{x) - u q{x) = j dyi dy2 P{x, 2/1,2/2) (108) 

dyi(l){x,yi,y2) dy^(t>{x,yi,y2) 
c{x) - 2u a{x) + u^q{x) (109) 

= j dyi dy2 P{x,yi,y2)[dy^(t>{x,yi,y2)f 

Q = {{m^)) (110) 
^ = 2Q-g(l) - - - + 

^ m-q[l)] ((^(tanh-^m-/?2/i))) (HI) 
C = Q + 2A + u{l-2Q + 2A)-u'^Q (112) 



where we have made use of Eq. (HO) in Eq. (HI) 
and Eqs. (110), (111) to yield Eq. (112). The average 
(((. . .))) is defined as 

((0(m,2/i,2/2))) = j dyi dy2 P(l,2/i,2/2) (H3) 
j^^ dm 0{m, 2/1,2/2) e-^f"" 



/Aa\' 



/^^ dm e'^'f"" 



or else, using Eq. (104), as 
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{{0{z,yuy2))) = j dyi dy^ P(l,yi,y2) (114) 

J^^VzO{z,yi,y2) p{-u,yi,z) e^^ 
I-oo'^^ P{-u,yi,z) 

The saddle point equations (110)-(112) for the elements 
of the diagonal matrices can, then, be written as 

g ^^tanh(/3 z v/A^ + /3 yi)^)) (115) 
Q Aa + A Aq = -Aq[l-Q{l+u)] (116) 

+ ^ ((^tanh(/J z^/A^ + p yi)] 

C = Q + 2A + u{l-2Q + 2A)-u^Q (117) 

_|,)_,(-;^ 

The Parisi Eq. (5) and the boundary condition Eq. 

(18) can be obtained as a BRST-SUSY rcdnction of Eqs. 
(98)-(99), setting A = a{x) = C = c{x) = 0, but leav- 
ing Q 7^ 9(1)- Furthermore, in our notation, q{l) be- 
comes q{xb) whereas Q is q{l\xb) as represented in the 
self-consistency Eq. (26). 

If a BRST-SUSY-breaking solution exists it must be a 
solution of the above equations with, as necessary (but 
not sufficient) features: a{x) ^ and/or c{x) 7^ 0, a lower 
band-edge equal to /eq and a positive xp decreasing to 
zero as / -» /eq. 

C. An instance of BRST-SUSY solution 

In Rcf. [23] the initial condition of the very same equa- 
tion (5) is given by Eq. (14), that we report here for 
clarity: 



Z{h 



(ii, 



v/27r[(9bm - r/(l)] \l-m 



, ^tanh ^ m — Am -|- hi^ 
i 2[Qbm - r?(l)] 

+ (^A - ^?7*(1)) m'^ + u f{m; 9) -h m /12} 



/(m;g) = - log2 - log(l - m^) 
1 

-|--mtanh"^ m - — (l - g|^) 



The above expression coincides with Eq. (80) if we set 
B ^0, Qbm = and A = A + p(l). The symbol / is 
what we call Pf in the present paper. 



Such a function is also identical to our boundary con- 
dition Eq. (18), provided a suitable change of notation 
is performed. First of all the a;-range has to be shifted as 



X € [0, 1] ^ X e [0, Xb] 



(118) 



so that any parameter computed in x = 1 will result 
as computed in Xb in our notation. Then we identify 
u = —Xb and on the line of transformations (90)-(92), 
(94)-(96) we make the following change of variables: 



Qbm = P''q{l\xb) 



V{x) = p\ix) (119) 



X = ^f3'q{l\xb) v*{x) = ^PMx) (120) 

A = ^p^q{l\xb) p{x) = -^p\{x){121) 

If Xb > Xc, q{l\xb) = Qea- Moreover the local fields are 
trasformed as in Eqs. (102)-(98), with yi = y and setting 
1/2 = without loss of generality: 



hi = -Py 



h2 = ^py (122) 



and the integration variable m (TAP-site-magnetization) 
is changed in z using Eq. (104). 

The function 0(x, yi, 2/2) involved is linearly connected 
to the solution of Eq. (5) with boundary condition (18) 
by the relation 

^BMy(x, y, 0) = /3X6<^Legendre(a;, 1/)H —{1-qif (123) 

Exploting Eqs. (118)-(122), (104) and (123), the com- 
plexity becomes equal to the Legendre transform of f''^^, 
given by Eq. (27). 

The actual computation performed by Bray, Moore 
and Young set, furthermore, the difference Aq = (Qbm — 
r]{l))/P'^ = q{l\xb) — qix'i^) equal to zero. 



D. Discussion 

The resolution of Eq. (98) keeping A, C, a{x) and 
c{x) different from zero from the beginning should lead 
to a BRST-SUSY-breaking complexity. Such a solu- 
tion would display some sclf-consistcntly determined or- 
der parameters, given by Eqs. (107)-(112), breaking the 
supersymmetry and allowing for the construction of a 
quenched complexity, function of some Xb break param- 
eter and its conjugated /. 

From what is already known about the annealed com- 
plexity and the property of marginal stability of the SK 
model at the lowest free energy values, wc expect this 
solution to fulfill precise requirements. First of all there 
cannot exist states of any kind below /eq, that is the 
minimal free energy value achievable by the system at 
a given temperature. Therefore the lower band edge /o. 
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that is less than /cq in the annealed approximation, must 
coincide with it in the quenched computation. Moreover 
the system is known to be marginally stable, i.e. the 

replicon, or Plcfka parameter, is zero at /cq. 

The BRST-SUSY-breaking annealed solution displays 
a strictly positive xp for any value of the /-conjugated u 
variable for which the complexity is positive. Numerical 
evaluations of this parameter, directly counting the num- 
ber of solutions of the TAP equations for a N spins sys- 
tem yield different, apparently contradictory, results. In 
Ref. [46] Bray and Moore performed the direct evaluation 
of the number of solutions of TAP equations and in Ref. 
[47] Plefka has carried out an analogous probe consider- 
ing modified TAP equations such that only solutions of 
the original equations with xp > arc selected. In both 
cases a finite size scaling analysis reveals that the repli- 
con eigenvalue/Plefka parameter goes to zero as N~^^^, 
thus indicating an exclusively marginal kind of stability. 
On the contrary, a very recent numerical study by Cav- 
agna et. at^^ displays an xp behaviour consistent with 
the one obtained for the annealed BRST-SUSY-breaking 
solution, hinting a strictly positive xp in the whole free 
energy support for which S > (sec Ref. [18] for the 
annelead complexity a;p). Even about the free energy 
density support numerical findings disagree, Plefka's re- 
sults hinting for a vanishing support^^ as N increases, 
whereas the results of Ref. [48] show a finite support 
(even consistent with the annealed support, at least for 
N < 80). 

Further analysis would be worthful to better under- 
stand the whole picture. In Ref. [9] the solution of the 
annealed saddle point equations was shown also to satisfy 
the quenched (replica symmetric) saddle point equations 
for all u above a certain, negative, value u. This includes 
u = 0, to which the largest value of the complexity cor- 
responds. Over the /-support this is equivalent to say 
that the annealed solution holds for / € [f{u),fth] but 
is unstable in the range [/o; /(«)]• Therefore, going to 
the exact quenched computation xp should decrease from 
xp{f{u)) > to xp(/eq) = in a well defined solution. 

The ensemble of TAP-solutions over which the saddle 
points are computed is of unknown nature, by this mean- 
ing that we cannot know a priori if the TAP stationarity 
points counted by our complexity are minima or saddles 
of any order and whether they lead to the correct linear 
susceptibility {xp > 0) or not (a;p < 0). As a conse- 
quence, we are not able to say whether we are counting all 
possible solutions (in which case a complexity breaking 
the supersymmetric invariance would be mathematical 
inconsistent) or only some solutions belonging to a par- 
ticular subset (in which case the invariance would not be 
a necessary requirement anymore). Recently Apelmeier 
et. have shown that the BRST-SUSY-breaking, else 
said BM,^ complexity is counting stationary points of the 
spin-glass free energy landscape, that are minima with 
one flat direction and that lead to xi = /^(l ~ q)- The 
flat direction is bounded to an isolated zero eigenvalue 
that can only be found considering subextensive correc- 



tions to the complexity saddle point. These very same 
corrections amoimt to a prcfactor of cxp NY,^^ that turns 
out to be zero at all orders in a 1 /N expansion when com- 
puted on the BRST-SUSY-breaking saddle point. This 
was shown by Kurchan"*^ for the case u = 0, correspond- 
ing to the maximum value of such complexity, (and to 
zero complexity in BRST-SUSY case). That argument 
can be extended to any value of u (or /), indicating that 
the prefactor can be either exactly zero or exponentially 
small in N.^^ Since the Morse theorem must hold even for 
the BRST-SUSY-breaking solution the vanishing prefac- 
tor should go like exp(— A^E^^) at u = 0. 

The zero prefactor and the isolated eigenvalue are dif- 
ferent ways of expressing the fact that metastable states 
in the SK model arc of a different nature with respect to 
those of the p-spin interacting model, with p > 3. An 
isolated eigenvalue has, indeed, been found also in the 
spherical p-spin model^ but there it is always positive. 
This difference is crucial to distinguish between the dy- 
namical behaviour of spin-glasses and structural glasses 
relaxation. In the case of p > 2, representing a good 
mean-field model for structural glasses, metastable states 
arc proper local minima (see e.g. [3 -5] for the spherical 
p-spin model or [20] for Ising p-spin) whereas in the SK 
model they turn out to be stationary points of different 
kinds in the thermodynamic limit. 



VII. CONCLUSIONS 

The main motivation for our work has been the need- 
ing for a clearer understanding of the two different ap- 
proaches that have been presented for the behavior of 
the complexity in the SK model up to now: the first one 
was originally introduced by Bray and Moore^ and the 
second one initially proposed by Parisi and Potters^^ and 
explicitely carried out (in a somehow different conceptual 
framework) quite recently in Ref. [14]. 

In Ref. [10] it was shown that the complexity could 
be obtained by calculating the partition function of 
TO distinct real replicas of the system^^ and provided 
the connection with the previous BM formalism by 
means of a generalization of the two-group An.satz}'^ In 
that framework they assumed an 'unbroken' two-group 
Ansatz in order to perform explicit computation. This 
has been finally performed exploiting the knowledge of 
BRST SUSY^'"^ (formally equivalent to the unbroken two- 
group) and leading to a supersymmetric complexity dif- 
ferent the BM one.i" 

The difference with the BM complexity has turned 
out not to be at the functional level, but simply the 
new solution is a second saddle point of the same 
functional. This means that there is just one BRST- 
SUSY functional, the BM complexity functional, dis- 
playing (at least) two saddle points, one satisfying the 
BRST-SUSY^* and the other one breaking it^. Both so- 
lutions identify an extensive complexity, computing the 
number of solutions of the Thouless-Anderson-Palmer 
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(TAP) equations^^ in the annealed approximation. In 
Ref. [18] the authors analyzed the two annealed solutions 
and their problems reaching the conclusion that both of 
them had to be discarded as candidate complexities of 
the SK model. 

In this paper we have faced the quenched calculation 
and the analysis of the properties of the complexity of the 
Sherrington-Kirkpatrick model making use of the Full 
Replica Symmetry Breaking Ansatz. Such complexity 
can be seen as the Legendre transform of the replicated 
free energy potential. A detailed description of the be- 
havior of the complexity is presented in Sec. IV. Re- 
sults are obtained both by numerical cxtremization of the 
replica free energy in the FRSB limit at a given temper- 
ature below the critical one (T^) and by anaytical expan- 
sion around the critical temperature. The complexity we 
find is the quenched analogue of the annealed complexity 
studied in Ref. [14]. As such, it satisfies BRST-SUSY in- 
variance. Even though going to the exact quenched com- 
putation cures some deficiency intrinsic in the annealed 
approximation (e.g. the lower band-edge now coincides 
with the equilibrium free energy) other problems arise 
or old inconsitcncies are not overcome. The most impor- 
tant outcome is that any attempt to build a BRST-SUSY 
complexity using the Legendre transform approach at fi- 
nite temperature leads to a negative replicon everywhere 
in the complexity support apart from the lower band- 
edge/equilibrium value /eq, where the replicon is zero 
and replica stability is marginal. This means that the 
solutions of the mean field TAP equations that our ob- 
servable is counting are not physically relevant. More 
specifically they do not express a local linear susceptibil- 
ity equal to /3(1 — qea)- This is related, in the replica 
language, to an instability against overlap fluctuations. 
Moreover, this leads to an even more serious problem, 
since a negative xp is mathematically inconsitent with 
the B = solution (see Sec. VI) of the saddle points 
equations for complexity, as shown in the Appendix of 
Ref. [18]. 

We are, then, brought to the conclusion that a com- 
plexity satisying BRST-SUSY, such as, e.g. the one of 



p-spin models {p > 2), cannot yield any extensive con- 
tribution in the SK model. A subextcnsive complexity 
was yield by the analysis performed by van Mourik and 
Coolen^^ in a dynamical reformulation of Parisi solution 
in terms of a multiple (infinite) hierarchy of decoupled 
time-scales, in which a small fraction of slower spins plays 
the role of an effective quenched disorder for the faster 
spins. Since in their construction the cjuantity of 'slow' 
spins turns out to be subextcnsive with respect to the 
quantity of 'fast' spins, the entropy of the slowest dy- 
namic variable, i.e. the comlexity, is zero. This means, 
though, that a "supersymmetric assumption" is hidden 
somewhere in their procedure. 

We eventually discussed the possibility of computing a 
BRST-SUSY-breaking quenched complexity, by means of 
the tools presented in Sec. VI. Even though the nature of 
the metastable states counted by such complexity would 
be different from the one of structural glass metastable 
states, and conceptually less intuitive, an extensive, ther- 
modynamically stable, complexity can, in principle, still 
be computed. 
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Appendix 

In this appendix we show the expansion of some ob- 
servables in the 'reduced' temperature r = Tg — T, where 
Tg = 1 is the critical temperature of the SK model. The 
following expansion are strongly asymptotics at higher 
orders, but, as already anticipated in Sec. IV B, Pade 
resummation allows for converging residts. 

The maximum complexity as a function of the reduced 
temperature is: 



2t^ 187 9938 2313541 10350722 r^^ 19122138382 r^^ 

j^max/ N ^ ^ 1 ^ /^24) 

^ ^ 81 81 1215 10935 382725 229635 51667875 ^ ^ 

57172333507 128525178819148 t" _ 182162380772635931 t^^ 133461584583918966554x^6 

17222625 ^ 3978426375 537087560625 34910691440625 

964957153958428975243 r^^ 591101435671536495956948 r^^ ^9 
20946414864375 999715254890625 -* 

The difference between the maximum value of quenched and annealed complexities is given by 

v^xnax. X _ vma.. ^ = ^ _ ^ , 3394 T^" _ 1381 t" 9052283 t^^ _ 788437 

^ ^ a w 243 135 54675 4725 7381125 273375 ^ ' 

6757091354 84666457147717 r^^ _ 136466739459370883x^*5 

217005075 107417512125 10741751212500 
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744513309305871199 50994683697175207152217 r^^ . 19. 
3989793307500 18851773377937500 ^'^ ^ 

The difference between the threshold free energies yield by quenched and annealed computation is 

f (^_fa(^_ t'^ 13r« 2711r9 70474t1" 1701148r" 453457547x^2 
hh{T) - /th('r) - " 810 ^ 54675 382725 ^ 5740875 ^ 82668600 ^ ' 

269738513263 779531487166381 27180847504347788 r^^ 

2411167500 477411165000 1220653546875 

211583950213407904679 13746420017310240661979 r^^ . i8\ 

698213828812500 3258331201125000 ^ 
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